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Abstract 

Presently, the most accurate value of the free neutron beta-decay lifetime result from measure- 
ments using fluid-coated ultacold neutron (UCN) storage bottles. The purpose of this work is to 
investigate the temperature dependent UCN loss rate from these storage systems. 

To verify that the surface properites of Fomblin films are the same as the bulk properties, we 
present experimental measurements of the properties of a liquid "Fomblin" surface obtained by the 
quasielastic scattering of laser light. The properties include the surface tension and viscosity as 
functions of temperature. The results are compared to measurements of the bulk fluid properties. 

We then calculate the upscattering rate of UCN from thermally excited surface capillary waves 
on the liquid surface and compare the results to experimental measurements of the UCN lifetime 
in Fomblin fluid- walled UCN storage bottles, and show that the excess storage loss rate for UCN 
energies near the Fomblin potential can be explained. The rapid temperature dependence of the 
Fomblin storage lifetime is explained by our analysis. 

PACS numbers: 14.20.Dh,61.12.-q,03.75-b 



I. INTRODUCTION 



The concept of a fluid-walled ultracold neutron (UCN) storage bottle is due to C. Bates 
who developed and demonstrated the technique at the Risley reactor in the U.K. The 
basic idea is that a hydrogen-free fluorinated oil (Fomblin) wets glass and metal surfaces, 
and thereby produces a clean, chemically perfect, microscopically fiat and reproducible UCN 
reflecting surface. Fomblin is quite viscous at room temperature; when it is sprayed onto 
surfaces, it drains away leaving a thin stable film. The fluid readily seals gaps and holes 
associated with UCN entrance and exit valves, provided that the mechanical parts fit rather 
closely. Fomblin also has a very low vapor pressure. Some of the most accurate determina- 
tions of the neutron /3— decay lifetime have been obtained with Fomblin-walled storage cells. 

There is some uncertainty in the chemical formula of Fomblin Y Vac 18/8 (manufactured 
by Ausimont/Montedison Group), the material most widely used in UCN storage experi- 
ments. The mean molecular weight is 2650, and the stoichiometry is roughly CsFgO, with 
density at 20° C of 1.89 g/cc, with a mean effective UCN potential of = 106.5 neV. The 
UCN loss coefficient (@], Eq. (2.66)) is given by the ratio of the imaginary W to real part 
V of the effective UCN potential, 

V Ana 

where cx; is the total 1/v loss cross section for neutrons with wavenumber k, and a is the 
coherent scattering length. / has contributions from nuclear absorption and bulk process, 
e.g. upscattering which occurs for temperature T > 0. For the range of formulae for Fomblin 
given in P, the possible values of / due to nuclear absorption lie in the range 

2.5 X 10'^ < / < 6 X 10"^ (1) 

The uncertainty in the Fomblin molecular formula given in [Q is a result of the coherent 
scattering lengths of C, F, and O being quite similar (6.65, 5.60, 5.81 fm, respectively) 
so the effective UCN potential V is relatively insensitive to the relative numbers of each 
atomic species. On the other hand, the nuclear absorption cross sections of C, F, and O are 
quite different (3.5, 9.6±.5, 0.19 mb, respectively, for thermal (2200 m/s) neutrons). Thus, 
the stoichiometry of Fomblin could be better determined by simultaneously considering / 
extrapolated to T = 0, which is quite sensitive to the atomic fraction of F because it has 
a comparatively larger absorption cross section, and V which cannot be used to effectively 
determine the relative amounts of F and O. 

For experiments operated near room temperature (290 K - 320 K), it is found that 

2 X 10"^ < / < 6 X 10"^ (2) 

and the temperature dependence of / is much faster than being proportional to T or ^/T as 
might be expected due to atomic and molecular fluctuations in the wall (see, e.g., 0], Sec. 
2.4.6) , and if fact seems to follow the viscosity temperature dependence rather closely. Some 
attempts were made to explain the observed temperature dependence of /. For example, 
if there is a 1/v cross section for slow neutrons due to upscattering in the fluid, this cross 
section can be used to modify W to include this loss mechanism and in this regard, the 
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transmission of 60 A neutrons through a sample of Fombhn as a function of sample temper- 
ature was measured [^. These measurements are discussed in where it is stated that the 
temperature dependence of / comes within a factor of 1.5 of explaining the experimentally 
observed temperature dependence of the loss rate. The loss rate for the lowest temperature 
218 K implies / = 6 ± 2 x 10^^ which is still an order of magnitude greater than Eq. 
(HI). More recent measurements by by Morozov et al. of the transmission of 9 m/s neutrons 
through a 4 mm thick sample as a function of temperature are consistent with these results 
i- 

It is possible that /(T) has contributions from different sources. The purpose of our work 
reported here is to evaluate the upscattering loss of UCN due to thermally excited surface 
capillary waves. Pokotilovski has investigated the energy broadening of (monochromatic) 
UCN stored in a fluid walled bottle |T0|; this is a question closely related to the upscat- 



tering loss, but it seems to us that the temperature dependent loss rate is an experimentally 
better studied phenomenon but is lacking a fundamental explanation. The principal data 
we will address were obtained by Richardson et al. |]TT . 



II. EXPERIMENTAL DETERMINATION OF THE THERMOMECHANICAL 
PROPERTIES OF FOMBLIN Y VAC 18/8 

Although the properties of this material are available from the manufacturer, we thought 
it prudent to perform measurements ourselves under conditions similar to those used in 
an experiment. The sample of Fomblin Y Vac 18/8 (hereafter referred to as Fomblin) was 
obtained from that used in These is anecdotal evidence that the viscosity changes 
for degassed samples in vacuum compared to samples that have been stored a long time 
in air. We therefore performed our measurements with the Fomblin under vacuum when 
possible, with the material being introduced into the various apparatuses to be discussed by 
distillation under high vacuum (P < 5 x 10~^ torr). 

These measurements were done in Aug. 1991; some details were lost (e.g., specific equip- 
ment models employed) but given the modest accuracy of the results, these details are not 
important. 



A. Mechanical Viscosity Measurement 

The viscosity was determined by the time it took for a known volume of Fomblin to drain 
through a Pyrex capillary of known length and diameter, under the influence of gravity. The 
glass system comprising two relatively large volumes connected by a capillary is shown in 
Fig. 1. Marks were placed on the large tubes corresponding to a volume v = 0.425 ± .005 
cm^. The capillary diameter was 2.00 it 0.01 mm, with length 4.8±0.1 cm. 

Fomblin was distilled into the glass measurement apparatus under high vacuum; the 
glass was sealed off from the vacuum system, and the time to drain the known volume was 
measured as a function of temperature, with the temperature determined by a water bath. 
The kinematic viscosity is given by (|1^, Eq. (17.10)) 



ng 



Kv/t) 



(3) 
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where v is the kinematic viscosity, g the acceleration of gravity, V is the known volume, t is 
the time to drain the known volume, and R is the capillary radius. The time to drain the 
volume varied from about 150 sec at the lowest temperatures, to about 2 sec at the highest 
temperatures; at the lowest temperatures, the constancy and value of the temperature was 
the major uncertainty (±1.2 K), while at the high temperatures, the timing (to 0.1 sec 
accuracy) was the principal limitation to the accuracy. In addition, there is a small correction 
associated with the change in height of the liquid layer as the Fomblin drains, but it is less 
than 5%. The results are shown in Fig. 2 and are in good agreement with, and follow the 
trends of, values reported by the manufacturer (1.9, .09, and .02 cm^/s at 293, 393, 493 K 
respectively). In the range of measurements, the viscosity is adequately described to within 
experimental error by the following empirical formula: 

v{T) = 2.0 X 10^6-° °^^^ + 900T-^/^ (4) 

where p is in cm^/s, and T is measured in K, with estimated accuracy over the temperature 
range 280-308 K which will be used later. 



B. Thermal Expansion 

Fomblin was distilled into a Pyrex capillary tube of 2 mm inside diameter, filling the tube 
to a length of about 5 cm. The tube was sealed off from the vacuum system, and the relative 
change in length of the sample as a function of temperature was measured. The results are 
shown in Fig. 3; the principal limits to accuracy were the length measurement (±0.2 mm) 
and temperature control (±0.5 K). The thermal expansion of Fomblin is about 100 times 
larger than Pyrex, so the expansion of the glass is an insignificant correction. The density 
as a function of temperature is given by 

p{T) - 1.89[1 ± 8.97 X 10-^(294 - T)] g/cm^ (5) 

where T is measured in K, and the accuracy of this fit is about 5%. The temperature 

dependence of the density is particularly important because it leads directly to a temperature 
dependence of the effective UCN potential. This information was not available from the 
manufacturer. 



C. Surface Tension 

The surface tension was measured using a Cenco student demonstration tensiometer 
which in essence measures the restoring force against a length of wire pulled up against the 
liquid surface. Unfortunately, there was no simple way to perform this measurement on a 
sample under vacuum. The apparatus was calibrated with deionized water. For Fomblin, 
the surface tension was found to be 

q; = 24 ± 1 dyne/cm (6) 

and was independent of temperature to within measurement accuracy. The value quoted by 
the manufacturer is 20 dyne/cm. 
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D. Liquid Properties by Light Scattering 



By use of a heterodyne laser scattering technique, the dynamic hquid surface properties 
of Fombhn could be measured directly without physically contacting the surface in any 
significant way. The idea is that the surface undergoes a sort of brownian motion due to 
thermal fluctuations. 

The theory describing the two-dimensional thermally excited surface capillary waves has 
been developed by Bouchiat and Meunier; when presented as a time and two-dimensional 
surface Fourier transforms, the mean-square height fluctuations are described by [pT^ 

1 



TXljJ p 



D{-iujT) 



(7) 



where 



D{S) = y + {1 + - V1T2S (8) 

y = ^^ r=-^ (9) 

where, as before, a is the surface tension, p is the density, and 77 is the viscosity rj = up 
where u is the kinematic viscosity. 

Directed laser light reflecting from this surface will be scattered by time-varying surface 
disturbances, with surface wavevector q and frequency u. By purposely aiming a fraction of 
the light specularly reflected from the surface in an angle corresponding to diffraction by a 
disturbance with wavenumber q (the "reference beam") to the detector, the time- dependent 
scattered component can be measured by beating with the reference beam. The reference 
beam is produced by placing a weak transmission diffraction grating in the specularly re- 
flected beam (this grating also intercepts the scattered light which is extremely weak). This 



technique was proposed and developed by Hard et al. \T^\ The specific apparatus that we 



used is described in |Ta], with the exception of the vacuum sample cell. As discussed in 



15| , data was aquired with the apparatus on a vibration isolated optical table. In addition, 
the vacuum sample cell isolated the liquid surface from disturbing air currents which were 
a problem for free liquid surfaces in air. 

The Fomblin sample was contained under vacuum in a cubical cell; the cell was con- 
structed of 5 mm thick float glass plates, and had internal dimension of 4 cm^, glued together 
with Torr-Seal Epoxy (Varian). A 5 mm diameter Pyrex tube was glued into a hole drilled 
in the top plate; Fomblin was distilled into the cell (to a depth of 3 mm) through this tube 
and then sealed off from the vacuum system. 

The cube was cooled/heated from the bottom with a Peltier device, and the temperature 
was controlled with a simple feedback loop. Light was directed onto the liquid surface at 
an angle of approximately 70° from normal through one of the vertical walls of the cell; the 
scattered and specularly reflected light exited the opposite vertical wall. The laser light was 
detected with a PIN photodiode, and the heterodyne signal was Fourier analyzed with a 
Hewlett-Packard audio frequency spectrum analyzer. Results as a function of temperature 
and q are shown in Fig. 4. Also shown are four-parameter fits to Eq. (|^); these parameters 
include r, y, an amplitude factor, and a DC offset. Results of the fits are given in Table I. 

The inferred values of u and a are also given in in Table I. As can be seen, these results 
agree with previous mechanical measurements; the values of u are plotted in Fig. 2, shown 
as squares. 



5 



III. CALCULATION OF THE UPSCATTERING LOSS RATE DUE TO 
THERMALLY-EXCITED SURFACE CAPILLARY WAVES 



The surface waves described in Sec. II D can inelastically scatter (diffract) a UCN that 
reflects from the surface. The probability to upscatter a UCN to an energy outside the 
UCN range can be calculated using Eq. (|^) when it is properly normalized. The formalism 
describing non-specular inelastic scattering was developed by Pokotilovski |TD[ and we will 
use his results with minor modifications. 



A. Normalization of Eq. (^) 

The work of [|T^ was intended to derive the spectrum of the surface fluctuations; the 
absolute normalization appears to have been of secondary importance. For our case, the 
overall magnitude is of crucial importance. 

First, the static average surface distortion, based on statistical mechanics considerations, 
is given as 



= ^ (10) 

where A is the surface area. If we imagine ^ to be a square with dimensions L so ^ = L^, 
the effective spread in and qy is 

Ag^ = Agj, = ^ ^ = dq^ dqy (11) 

and therefore, to cast Eq. (^ per unit qxqy, a factor of l/(27r)^ is required. 
Second, the integral of P{u,q) over lj should give the static result, Eq. (p!0|): 



P{uj,q)duj = — — -r^Im 

TT p 



dx 


1 


/ — CO 
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r - = (12) 

IT p y aq^^ 



where the integration was performed by noting that the only pole of the integrand is at 
X = 0; by integrating along the real axis and taking a semicircle around the origin, the 
integral is —m times the residue at x = which is l/{—iy). 



B. Possible Corrections to Eq. (0) 

The coating obtained by spraying Fomblin and allowing the material to drain to a thin, 
stable film are of order 10"*^ cm thick or greater. The multiplicative correction factor to the 
capillary wave dispersion relationship, due to finite film thickness, is very roughly propor- 
tional to tanh(g(i) where d is the film thickness (|l^. Sec. 62, Problem 1). For q of interest 



for UCN scattering, 10~^cm x q » 1 so tanhgc? ^ 1 so we expect no significant effect due to 
film thickness. However, at higher temperatures where the viscosity is low, the films might 
drain to layers thin enough so that there is a substantial correction. The net effect will be 
to reduce the upscatter rate. 
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We are also considering rather high frequency waves; it might be expected that the 
Fombhn has some frequency dependence. However, Fombhn appears to be a Newtonian 
fluid (e.g., light scattering and mechanical measurements had no anomalous behavior) and 
the Fomblin molecules are not very large on the scale of what is considered a polymer. The 
molecular collisional frequency of Fomblin molecules, which is on of the factor that determine 
the viscosity, is much higher than the frequencies of interest for UCN scattering (e.g., 10^ 
Hz compared to more than 10^° Hz). For the analysis presented here, we will assume that 
the low frequency measurements described earlier in this paper are applicable. 



C. Kinematics 

Let us consider an incident neutron with wavevector (we set the azimuthal angle of the 
incoming wave to zero without loss of generality); referring to Fig. 5, 

ki = kixX + kizZ — ki{sm9iX + cos 9i z); — 2mEi/h^ (13) 

where the polar angle 9i is defined relative to the surface normal which is taken along z, and 
Ei is the incident UCN energy. The surface wave disturbance can be described by 

q^x + Qyy ^ q{cos(f)gX + sin 4>gy). (14) 

where is the azimuthal angle relative to the incident neutron momentum. 

When a UCN diffracts from the time varying surface, energy and momentum are con- 
served (we consider first upscattering, e.g., the UCN gains energy). 

= ^ + (15) 

We can consider a UCN of specified ki, hence |ki| is specified, and choose u which 
determines |kf|; choosing q determines Of and 0j, the polar angles of the outgoing neutron 
wavevector. The components of the wavevector in the liquid surface plane (denoted by 
subscript s) change on reflection by 

kfs = kis + q (16) 

or 

kfx — ki sin. 6 q cos (f)] kfy — qsincf). (17) 
For the component of kf normal to the surface, 

= ^/ - ^% - ^ (18) 

where the inequality constrains the values of q for allowed scattering. This sets a range on 
the magnitude of q: 

< g < qmax ^kf + ki. (19) 
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In particular, we can define a kinematic factor sucli tliat 



K{Ei,9i,uj,q, 



1 if k}^ > 
otlierwise, 



(20) 



wliere kj^ is given by Eq. 



D. Loss Probability per Bounce 



Tlie probability of loss per bounce was calculated by Pokotilovski |1T0| : 

kiz iKiz 



w. 



g,UJq 



8kizkf^P{uj,q) 



kfz + iK 



where 



8Fkizkf^P{uj,q) 



(21) 



(22) 



with kc the critical momentum for UCN reflection. This equation results from determining 
the scattered wavefunction amplitude then calculating the probability current normal to the 
surface. Pokotilovski was interested in below-barrier upscattering in which case the factor 
F = l. 

In the case where the final UCN energy is higher than the wall potential, some care must 
be used in calculating the upscatter loss probability. The UCN will be lost from the system 
if the final neutron energy is greater than the wall potential. This loss has two components: 
1. The upscattered reflected wave will be lost from the system, 2. The component of the 
scattered wave within the wall (film) is not an evanescent wave, but is a propagating wave 
(e.g., Kfz in Eq. (6) of [|TU| or Eq. (^) above is imaginary). We must therefore account 
for this losses for the wave transmitted into the material. Referring to Eqs. (A6.21) and 
(A6.22) of [Q, the amplitudes of the reflected and transmitted upscattered waves are equal. 

in 



The probability current perpendicular to the liquid surface determines the loss 0], e.g. 
7l| the factor kfz ^ kfz + \ Ikl^ — /c?. Thus, 



Eq. 



if k% < kl 
otherwise. 



(23) 



The UCN will be lost from the system if its final energy is greated than the wall potential 
V . We can recast this in terms of allowed wave frequencies: For UCN loss to occur, 

uj + ujq> uJc = Vlh; fkjQ = Ei (24) 
where Ei is the incident UCN energy. Thus, there is a minimum uj for upscattering loss. 



(25) 
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The differential probability that a UCN, with Ei, 6i, and oo specified, is upscattered and 
lost from the system due to surface waves in a small range around a value q is given by 

dfisw = SFki^kf^Wi ^^^^^ dq^dqyduj (26) 

where Wi is a weighting factor giving the probability for a UCN to be incident at angle 9i, 

Wi = 2 sin 9i cos 9i, (27) 

and dqxdqy = qdqdcj). 

The total probability for loss is equal to the integral of dfi over q and (j), followed by an 
integration over 6i, and finally over uj with uj > Umin, 



IJswiEi,T,UJmin) = / duj WidOi 2 d(f) 

J UJrnin JO L ^0 



P(uj a)' 

qdq K{Ei, 9i, u, q, (fySFki^kf^ \ 

[2txY _ 

(28) 

where K is defined in Eq. (pO]), and the dependence of iisw on the temperature (due to, e.g., 
the viscosity and density change of Fomblin) and incident UCN energy Ei and minimum 
energy uJmin to upscatter to greater than V are explicitly indicated. 

The integration is analytically unwieldy, so it has been performed numerically. The 
results, taking into account the temperature dependence of the viscosity and density (hence 
effective UCN potential) of Fomblin are shown in Fig. 6. 

The slow curvature with increasing temperature (hence decreasing viscosity) results from 
the y term in Eq. (|^, which is proportional to the surface tension, becoming important. 



E. Reanalysis of Previous Data 



It is noted in [11| that a simple fit to a value of / and to losses through holes and gaps 
of total area Ah does not explain the observed lifetimes of a fluid-walled bottle for energies 
within 10 neV of V . There is enough flexibility in / and to adequately flt the lower 
values of energy; as seen in Fig. 6, the surface wave upscatter loss probability drops rapidly 
with energy. However, the temperature variation in / as derived from these flts seems 
non-physical. 

In general, a full monte-carlo analysis of the system is required. For example, in a realistic 
storage system, the weight function Eq. ( P?] ) becomes a function of geometry. However, 80% 
of the UCN loss for Eq. (|28|) occurs for angles between 20° and 60°, so at the 10-20% level 
of accuracy, we can neglect modiflcation of the weight function. Furthermore, a full monte- 
carlo analysis would require knowledge of the liquid surface specularity for UCN reflection 
which is not known. Another fortunate feature of the Ref. measurements is that the 



cylindrical storage cell height 2R = 6.6 cm leads to a small change in UCN energy as a 
function of position in the cell (7% effect) that can be modelled, as described below, with 
10% accuracy. 

For our reanalysis of the Ref. |]ll| experimental data, we use oure experimentally deter- 
mined values for viscosity (Eq. (4)), density (Eq. (5)), and surface tension (Eq. (6)), and 
neglect the uncertainties in these values. We fltted the observed lifetime vs. energy data to 
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/ and Ah was done as in [|TT|, but included all data points and the surface wave upscatter 
loss described by Eq. (]28|). However, particularly for Ei close to V, we need to take into 
account the effects of gravity in the storage bottle which was a horizontal cylinder with 
radius R = 3.3 cm and length L = 160). There are two effects that need to be accounted 
for. First, the upscatter loss function energy dependence is very steep for E^ approching V. 
Second, a UCN that upscatters at the "top" of the storage vessel will gain energy as it falls, 
and will eventually be lost from the system if the final energy Ej > V — 2mgR |]TB|. This 
effect can be modelled by varying Umin in Eq. ( P^D as a function of height in the bottle. 

The total average lifetime for a UCN with initial energy E (monochromatic) specified 
at a horizontal plane that intersects and divides the storage vessel into two halves at the 
cylinder axis can be readily determined. The rate of loss can be written 



where rg is the neutron beta decay lifetime, £ is the mean free path, 

2L 2R 



(29) 



(30) 



v{E) is magnitude of the UCN velocity at the plane, and fitot{E) is the total loss per bounce 
averaged over the storage vessel walls. The height- dependent total loss is given by 



l^tot{E, 



mgR cos 9 
E 



1/2 



Ah ' 

n{E — mgRcosO) + fJ^swiE — mgR cos 9, T,LJrnin — 2mgRcos9) + — — 

4:Vb_ 

(31) 



where 9 is the cylinder polar angle with 9 = being the bottom of the storage vessel, and 
/x(£') is the 1/v material loss term (0], Eq. (2.70)), 



ft{E) = 2f 



V 
E 



arcsm 



1/2 



V 
E 



1/2 



(32) 



As described in the Introduction, / has contributions from both nuclear absorption and 
upscattering and represents the losses from the UCN evanescent wave within the film during 
the time of reflection. 



The first factor in Eq. (^) expresses the change in wall collision frequency as a function 
of height and has a 2% effect for E ^ V. The average of the storage walls is given by 



l^tot{E) 



vr 



fitot{E,9)d9. 



(33) 



This average was numerically computed and then averaged over with the the energy resolu- 
tion function given in |Tl|, Eq. (10), to produce total loss curves as a function of temperature 
T; / and Ah are left as fit parameters for a given T. 

The most extensive data presented in [|1^ is for the temperature T = 294 K. Including 
the surface wave upscattering loss and fitting to the two parameters results in an excellent 
fit, as shown in Fig. 7. The decrease in (see Table II) compared with the fit results 
presented in [|lT| shows the importance of the surface wave loss for UCN energies near the 
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wall potential. It should be noted that the excess loss is not associated with slow heating of 
the UCN, as suggested in but with a direct upscattering loss. As will be shown in the 
next section, the heating/cooling effect probably does not alter the functional form shown 
in Fig. 7. 

Fits to the 283 K and 308 K data are shown in Figs. 8 and 9, with fit results tabulated 
in Table II. Again, including the surface wave upscattering loss provides an explanation of 
the extra loss rate for the higher energy UCN. 

It should be noted that the for these fits is extremely sensitive to the energy, e.g., 
changing the energy for either of the two highest points in Fig. 7 by ±1 neV changes 
by a factor of 0.5 to 2. Therefore, the fact that the spectrum shape can be varying during 
storage might be evident in the data and analysis, but as can be seen by directly considering 
Figs. 7-9, the simple analysis presented here is adequate to within statistics, implying an 
accuracy of 10-20%. Including the uncertainties in the measured Fomblin properties would 
lead to a modest reduction in x^. 

Finally, the values of / for the various temperatures shown in Table II indicate a signifi- 
cant temperature dependence. These values, along with the value derived from discussed 
in the Introduction are plotted in Fig. 10. It can be seen that a crude extrapolation to 
T = implies a reasonable range of values for the residual (nuclear absorption) loss. 

The observed variation of bottle lifetime as a function of temperature is about 4 s/K 0; 



referring to Fig. 7, the contribution to the loss rate from Eq. ( |32D is about a factor of two 
larger than the surface wave loss rate for 70 neV UCN. Interpolating from Fig. 6, the change 
in is about 5 x 10~^/K for 70 neV with the Fomblin near room temperature, while from 
Table II / varies as df /dT = 2.3 ± 1 x 10~^/K. The variation in r with change loss rate(s) 
can be estimated from Eqs. (27) and (30) (the effect due to the change in the potential V 
from thermal expansion is negligible) 



6t/5T = T^6fitot/6T = T^j [1.445/ + 5/^,^] /5T (34) 

which, for a 20 cm mean free path for 70 neV neutrons, with a storage lifetime of 710 s 
implies 

5t/5T = (710 s)'(1.44 x (2.3 ± 1) x 10"' + 5 x IQ-^)^!^^^ / ST = 3.6 ± 1.6 s/K (35) 

which is in good agreement with 4 s/K. The surface wave loss amounts to 20% of the net 
total loss for 70 neV UCN at T = 294 K. In [Q, it is stated that the observed upscatter loss 
in the transmission of 60 A neutrons comes within a factor of 1.5 of explaining the observed 
loss. Within the approximations here, it would seem that the surface wave scattering can 
account for most of the excess loss observed in this storage experiment. 



F. UCN Heating 

Heating of UCN occurs when energy is gained as before, but when the final energy is still 
lower than the wall potential. The calculation is carried out by integrating from near u = 
to LOmin as defined in relation to Eq. (pSf ). The lower limit is set by the resolution of the 
measurement apparatus, or, e.g., by the spectral width of the stored nearly monochromatic 
UCN, taken here as ±2 neV, so ii^^ = 2 neV and Ur = Er/h. The heating probability is 
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thus obtained by modifying Eq. ( ]28|) (e.g., the integration over oj is not performed to get 

PheatiEi,T,UJr) = llsw{,Ei,T , uj)duj (36) 

J UJr 

where uJmin was defined in relation to Eq. (|28|) and in this case represents the maximum 
energy change that a UCN can receive and remain trapped. The average energy change is 
given by 

pUJrnin 

AE{Ei,T,UJr) = tu^llsw{Ei,T,u)duj/Pheat{.Ei,T,Ur). (37) 

J UJr 



It is suggested in |TT| that energy might be transferred to the UCN at a rate of order 
SErms = 2(1.5) X 10~^^ eV per colhsion. The numerical results shown in Fig. 11 imply an 
energy change of around 10^^'^ eV per collision, obtained by taking the upscatter probability 
times the energy change. Of course, very small energy changes (less than the instrumental 
resolution) are more probable, but many more reflections are needed to move a UCN from 
the initial energy distribution; the 10~^^ eV/coUision provides a reasonable estimate of the 
heating rate. This result is in agreement with the results of Pokotilovski ^ and we thus 
do not further elaborate the heating issue which is more fully discussed by him. 

More recent measurements by Bondarenko et al. fl^ indicate a heating probabilities 
shown in Table III. In this experiment, UCN in a 0-52 neV energy range were stored in a 
Fomblin-coated bottle, and upscattering to energies between 52 and 106 neV were detected. 
To calculate this probability, we assume an initial UCN energies in the 0-52 neV range, 
and calculate the total probabilities to scatter to the 52-106 neV range, using Eq. (p6|). 
The energy-dependent upscatter probabilities are then averaged over the initial spectrum 



presented in |jT^, Fig. 5. As presented in Table III, it can be seen that our results are in 



reasonable agreement with the experimental results, to within the reported factor of 2 to 3 
uncertainty associated with neutron transport efficiency in the experimental apparatus. A 
plot of the upscatter probability as a function of temperature, for the parameters of this 
experiment, is shown in Fig. 12. 

The cooling rate and average energy decrease can be calculated in a similar fashion. In 
this case, the integration limits of Eq. ( p6| ) are and Ei/h — Ur- The probability to upscatter 
or downscatter is shown in Fig. 13. As discussed by Pokotilovski, the functions are very 
symmetric about the incident UCN energy, so the effect of a small energy change, either an 
increase or a decrease, tend to be averaged away. 

The measurements by Bondarenko et al. also determined the probability for UCN 



with initial energy around 13 neV, to downscatter to final energy in the range 0-11.5 neV, 
and found a probability of about 10~^ per bounce. The temperature of the Fomblin was not 
specified; however, the measured probability is consistent with the results presented in Fig. 
13 where the integrated probability to downscatter from 10 neV initial energy to less than 
8.5 neV is in the 10~^ range for both temperature curves shown in Fig. 13. (The variation 
in integrated probability with initial energy is fairly slow; the calculated result would not 
be significantly different for 13 neV initial energy.) The rapid drop in probablity with 
decreasing maximum final energy that is evident in Fig. 13 is consistent, within statisitics 



and uncertainties in neutron transport, with the data presented in |17|, Fig. 2. 
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IV. CONCLUSION 



As shown in Figs. 7-9, the UCN hfetime for energies near 100 neV (e.g., the high energy 
data points presented in Ref. |jll|) are well- described by upscattering loss due to surface 
waves. In this analysis, the low-frequency kinematic viscosity and surface tension were used. 
The excellent agreement between theory and experiment suggest that no modification of the 
low frequency dispersion curve for the surface waves are required, but we might anticipate 
that for higher accuracy data, this simple analysis ultimately might not be sufficient. 

It was suggested in []1T| that the excess high energy loss is due to changes in the stored 
UCN spectral properties; the analysis presented here indicates that direct upscattering to 
energies higher than the wall potential fully account for the observed loss rate beyond that 
due to nuclear absorption and upscattering in the film (as parameterized by / in Eq. (P^)) 
and losses through holes and gaps. The decrease in foi' the re-analysis of the Richardson 
et al. data, with no adjustable parameters, indicates that surface wave upscattering laregely 
describes the excess loss rate for UCN energy close to the wall potential. As described in 
Sec. F, and shown Figs. 11 and 12, the heating and cooling probabilities are sufficiently 
small so that, at the level of accuracy of the data presented in |]ll|, the effects of spectral 
evolution are not important. However, such effect are likely important for high- accuracy 
determination of the neutron lifetime from Fomblin coated bottles. 

We have shown that the observed temperature dependence of the net loss rate has a 
significant component due to the surface wave upscattering loss process. In , the measured 
temperature dependence of the lifetime for 70 neV in a Fomblin coated bottle is given as 
4 s/K. This was not fully explained by the temperature dependent 1/v cross section (e.g., 
the temperature dependence of / in Eq. (|32D). The temperature dependent loss implied by 
measurements of transmission of 60 A neutrons through Fomblin was 33% too small 
However, from our analysis, the surface wave upscattering accounts for up to 20% of the loss 
for 70 neV UCN. This accounts for much of the reported discrepancy. This, together with 
our estimate of the temperature variation of /, explains the UCN loss in Fomblin coated 
bottles to within experimental errors. 

Recent Fomblin transmission measurements using 9 m/s neutrons by Morozov et al. 
1^ indicate that df /dT = 9 ± 5 x 10^^/K which is in rough agreement with the result 
presented in Sec. HE, based on Fig. 6. Furthermore, the values of / implied by the 
Ref. P] measurements are in agreement with the results of our analyis presented in column 
five of Table II. Unfortuantely, our value of / extrapolated to T = (Fig. 10) has too 
much uncertainty to help with the analysis of the Fomblin stoichiometry as discussed in the 
Introduction. 

We hope that the analysis presented here will stimulate high-accuracy experimental mea- 
surement of the temperature dependent 1/v cross section of Fomblin. It is our expectation 
that the temperature dependence of / together with the surface wave loss rate can fully 
explain the temperature dependence of the storage lifetime of Fomblin coated bottles. 

We thank Albert Steyerl for critically reading the manuscript and providing clarifications 
on several points. SKL was supported by LANL |ldrd-dr 2001526. 
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TABLE I: Results of fits to the data shown in Fig. 4. 
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TABLE II: Comparison of UCN lifetime fit results without surface wave loss (Richardson et al, 
[9]) to results using the surface wave upscattering probability given in Sec. HIE. In both cases, 
the number of fit parameters are the same (two, / and A^). The results from [9] do not include 
the higher-energy lifetime data points; given in [9] has been adjusted here to include all points 
(those near 100 neV were excluded from the fit in [9]). 



Ref. [9] fit results 


With Surface Waves 


Temp [K] 


/ [10-5] 


Ah [mm^] 


(d.o.f.) 


/ [10-5] 


Ah [mm^] 


X^ (d.o.f.) 


283 


1.9(0.2) 


3.0(0.5) 


5(4) 


0.83(0.10) 


5.6(0.5) 


1.5(5) 


294 


2.45(0.1) 


0.0(0.0) 


8(8) 


1.28(0.10) 


3.5(0.6) 


2.1(8) 


308 


3.9(0.5) 


3.4(1.5) 


9(2) 


1.42(0.3) 


8.3(1.0) 


2.4(3) 



TABLE IIL Recent experimental measurements of UCN upscatter probability [15] compared with 
calculation. 



Temp (K) 


Exp. Upscatter Prob. 


Calc. Prob. 


343 


9 ± 1 X 10-6 


2.4 X 10-5 


308 


3 ± .8 X 10-6 


5.8 X 10-6 


298 


1.2 ± .8 X 10-6 


3.4 X 10-6 



Fomblin 



AVacuum Tip-Off 



Fiducial 
Marks 



48 mm 



2 mm 



FIG. 1: Schematic of the viscosity measurement apparatus. 
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Empirical fit: v(T)=2x10V° °^^"^ + 900/T^^^ cm^/s 
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FIG. 2: Temperature dependence of Fomblin viscosity from mechanical measurements (crosses) 

and laser light scattering (squares). 
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FIG. 3: Volume expansion of Fomblin with temperature. 
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FIG. 4: Experimental measurements of heterodyne detected inelastic light scattering with theoret- 
ical fits, for various temperatures as indicated in plots. 




FIG. 5: Kinematics of UGN scattering from surface waves; final energy greater than initial energy. 
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FIG. 6: Probability for UCN to upscatter to energy higher than the wall potential, as a function of 

temperature. Included arc the temperature dependencies of the viscosity and of the density. The 
effective potential is modified for the density change. 
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FIG. 7: Least square fit to 294 K data from [11]. The experimentally determined values p = 1.89 
g/cc, V = 2.0cm^/s, and a = 24 dyne/cm^ were used in the fit to and /, resulting in a reduced 
= 2.1 for the fit. Upper plot shows the net lifetime, while the lower plot shows the contributions 
from the various loss mechanisms. 
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FIG. 8: Least square fit to 308 K data from [11]. The experimentally determined values p = 1.87 
g/cc, V = l.Ocm^/s, and a = 24 dyne/cm^ were used in the fit to and /, resulting in a reduced 
= 2.4 for the fit. Upper plot shows the net lifetime, while the lower plot shows the contributions 
to the loss rate from the various mechanisms. 
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FIG. 9: Least square fit to 283 K data from [11]. The experimentally determined values p = 1.91 
g/cc, V = 3.0 cm^/s, and a = 24 dyne/cm^= were used in the fit to and /, resulting in a 
reduced = 1-5 for the fit. Upper plot shows the net lifetime, while the lower plot shows the 
contribution to the loss rate from the various mechanisms. 
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FIG. 10: Loss factor from the surface wave least-square fits (column five in Table II) and the 218 
K value derived from [7] plotted as a function of temperature. A simple extrapolation to T = 
yields a range of / consistent with the nuclear absorption value. 
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FIG. 11: Upscatter heating probability, with final energy 2 neV greater than the initial energy. 
Lower plot, the average change in energy with upscattering. 



24 



X 10" 




240 250 260 270 280 290 300 310 

Fomblin Temperature [K] 



320 330 340 




280 290 300 310 
Fomblin Temperature [K] 



FIG. 12: Upper Plot: The upscatter probability per bounce for UCN in the energy range 0-52 neV 
range, with initial distribution given in Fig. 5 of [17], to upscatter to the 52-106 neV energy range, 
as a function of Fomblin temperature. Lower Plot: The upscatter probability times the viscosity, 
divided by temperature, which is approximately constant for T < 300 K. This is because in this 
temperature range the y term in Eq. (8) is relatively unimportant, and S « 1. This leaves Eq. 
(8) proportional to T/r){T), with the integration over the q and uj factors contributing an overall 
factor roughly independent of T. 
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FIG. 13: The probability for UCN to up- or down-scatter, per neV final energy, for three different 
UCN initial energies. These results are similar to those presented in [9,10]. 
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